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ABSTRACT 

Advanced LIGO and Advanced Virgo will be all-sky monitors for merging compact objects within a few hundred 
megaparsecs. Finding the electromagnetic counterparts to these events will require an understanding of the transient 
sky at low redshift (z < 0.1). We performed a systematic search for extragalactic, low redshift, transient events in 
the XMM-Newton Slew Survey. In a flux limited sample, we found that highly variable objects comprised 10% of 
the sample, and that of these, 10% were spatially coincident with cataloged optical galaxies. This led to 4 x 1 0 4 
transients per square degree above a flux threshold of 3 x 10 -12 erg cm -2 s -1 (0.2-2 keV) which might be confused 
with LIGO/Virgo counterparts. This represents the first extragalactic measurement of the soft X-ray transient 
rate within the Advanced LIGO/Virgo horizon. Our search revealed six objects that were spatially coincident 
with previously cataloged galaxies, lacked evidence for optical active galactic nuclei, displayed high luminosities 
~10 43 erg s _1 , and varied in flux by more than a factor of 10 when compared with the ROSAT All-Sky Survey. At 
least four of these displayed properties consistent with previously observed tidal disruption events. 
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1. INTRODUCTION 

The X-ray band provides an opportunity to find high confi- 
dence counterparts to the compact object mergers that will be 
discovered with the second generation LIGO and Virgo grav- 
itational wave detectors. Within this decade Advanced LIGO 
(Harry et al. 2010) and Advanced Virgo' 1 (Acernese et al. 2008) 
are expected to begin detecting mergers of binary neutron stars 
and neutron stars with stellar mass black holes out to distances 
of a few hundred megaparsecs (Abadie et al. 2010). Placing 
these mergers in an astrophysical context and maximizing the 
scientific returns will require finding electromagnetic counter- 
parts to the events (Bloom et al. 2009; Phinney 2009). However, 
the positional accuracy of the gravitational wave detectors will 
be limited to tens or hundreds of square degrees (Abadie et al. 
2012; Fairhurst 2009; Klimenko et al. 2011; Nissanke et al. 
2013). Thus, associating an electromagnetic counterpart with 
a LIGO/Virgo detection will require an understanding that a 
chance coincidence within the LIGO/Virgo horizon is unlikely, 
even within a large sky region. In the optical band, large area 
survey instruments such as Pan- STARRS, Palomar Transient 
Factory, SkyMapper, and the future LSST will have a daunting 
challenge separating LIGO counterparts from stellar variability, 
supernovae, and other confusion sources (Kulkarni & Kasliwal 
2009; Nissanke et al. 2013), but may leverage the large ongoing 
effort to create schemes of automated or semi-automated tran- 
sient classification (e.g., Bloom et al. 2012). The radio band may 
also be searched for counterparts to GW events (Predoi et al. 
2010; Lazio et al. 2012), and large area searches for radio tran- 
sients are rapidly developing (Stappers et al. 2011; Bhat et al. 
2013). 

In the X-ray band, low number counts at bright flux levels 
may make identification of a LIGO/Virgo counterpart more 


6 https://tds.ego-gw.it/itf/tds/file.php ?callFile= VIR-0027A-09.pdf 


straightforward. However, studies of X-ray variability have 
tended to focus on persistent or repeating sources, particularly 
active galactic nucleus (AGN), X-ray binaries, or stellar flares. 
Past all-sky searches for soft X-ray transients include searches 
in the ROSAT All-Sky Survey (RASS) data (Greiner et al. 2000a; 
Fuhrmeister & Schmitt 2003), work with the XMM-Newton Slew 
Survey (Esquej et al. 2007; Starling et al. 2011), and searches 
for flashes lasting a few seconds (Gotthelf et al. 1996; Connors 
et al. 1986). 

A small fraction of compact object mergers are thought to 
create short gamma-ray bursts (GRBs; e.g.. Fox et al. 2005; 
Eichler et al. 1989), and if one occurred within the Advanced 
LIGO/Virgo horizon, the X-ray afterglow could be bright 
(~ 10 12 erg cm -2 s _1 ) for about a day after the merger (Kanner 
et al. 2012). Short GRBs may be associated with mergers of 
two neutron stars, or a neutron star with a black hole, and are 
typically characterized by a prompt emission lasting less than 
2 s, and a spectrum that is somewhat harder than the more 
prevalent “long GRBs” associated with stellar core collapse 
(Nakar 2007). Some double neutron star mergers may also 
exhibit a bright X-ray counterpart due to emission from a 
magnetar-powered ejecta (Gao et al. 2013; Zhang 2013; Metzger 
et al. 2008), observable for up to a few thousand seconds after 
the merger. Such signals could be detectable with a wide field 
focusing instrument such as the proposed ISS-Lobster (Camp 
et al. 2013) or A-STAR (Osborne et al. 2013), or in some 
cases with multiple observations of Swift /X RT (Kanner et al. 
2012; Evans et al. 2012). In order to estimate the ability to 
identify a unique source in such a campaign, we have searched 
the XMM-Newton Slew Survey Clean Source Catalog, version 
1.5 (XMMSL1; Saxton et al. (2008)) for objects consistent 
with counterparts to compact object mergers observable with 
the future Advanced LIGO and Advanced Virgo observatories. 
We sought both to measure the density of such events that are 
observable in a given moment and to characterize the nature 
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Figure 1. The fraction of XMMSL1 objects listed with RASS matches in each 
flux bin. The vertical line indicates the flux limit selected for the search. 


of the objects we found. Our search criteria emphasized low 
redshift, transient objects, and were similar to those of Esquej 
et al. (2007). However, this study used a data set covering five 
times the slew survey area available in 2007 and included a 
systematic measurement of the transient density using a simple, 
easy to emulate definition. 

2. CANDIDATE SELECTION 

We designed our selection criteria to seek objects consistent 
with a transient event within the Advanced LIGO/ Virgo neutron 
star merger horizon distance. We used the XMMSL1 soft band 
of 0.2-2 keV for all flux measurements, since it is similar to 
the ROSAT PSPC band. The XMMSL1 includes identifications 
with RASS sources with a 30"search radius (Voges et al. 1999; 
Saxton et al. 2008). We found that 80% of XMMSL1 sources 
brighter than 3 x 1 0” 1 2 erg s 1 cm -2 had matches in the RASS 
with no improvement in overlap for brighter sources, so we 
took this as the flux limit for our search (see Figure 1). For each 
XMMSL1 object brighter than this threshold, we attempted to 
place a flux upper-limit in the RASS data set. We explored a 
range of different radii for the source extraction region using 
data corresponding to ROSAT detected sources, and found that 
a 205" radius was needed to recover the median source with 
90% of the expected flux. For the XMMSL1 objects above our 
flux threshold, Bayesian upper-limits corresponding to a 2 a 
confidence level were applied to the ROSAT PSPC counts found 
within a 205" extraction radius. To convert from counts to flux in 
the 0.2-2 keV band, we used webPIMMS 7 8 with the same source 
assumptions that were used in the XMMSL1, namely, a power- 
law spectra with an index of 1 .7 and a hydrogen column density 
of 3 x 10 20 cm” 2 . For each object, we took a ratio between 
XMMSL1 flux as listed in the catalog and the RASS upper-limit 
we constructed, and kept only objects with a flux ratio greater 
than 10 (see Figure 2). These objects, or any objects with a 
flux ratio greater than 10, are referred to as “transients” for the 
rest of this paper. Since we were interested only in extragalactic 
objects within a few hundred megaparsecs, we expected their 
host galaxies to be visible in large optical and near infrared 
surveys such as the Sloan Digital Sky Survey s (SDSS) and 
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Figure 2. Sources above our flux threshold in the XMMSL1. The y-axis 
represents the flux ratio between the XMMSL1 and RASS. All RASS fluxes are 
calculated as 2<r upper limits, with red triangles indicating a number of RASS 
counts consistent with a non-detection at the 2 a level, and blue circles indicating 
a source was detected. A yellow, hollow circle indicates an identification with a 
galaxy in the XMMSL1. 

(A color version of this figure is available in the online journal.) 

the Two Micron All Sky Survey. Therefore, we demanded that 
candidates be listed in XMMSL1 as spatially coincident with a 
known galaxy or galaxy cluster with a 30" match radius. Finally, 
we checked our candidates for an AGN association, and for any 
other observations in the HEASARC database. 9 

3. RESULTS 

3.1. A Snapshot of the Sky 

Our search showed that source variability seemed to naturally 
divide the XMMSL1 survey into two classes (see Figure 2). The 
first class, representing roughly 90% of sources, was mainly 
below the flux ratio threshold value of 10, and mainly detected 
in both the RASS and the XMMSL1 , indicated by blue circles in 
the figure. The distribution of sources detected in both XMMSL 1 
and the RASS is shown in Figure 3, and is loosely consistent with 
a lognormal distribution. Their flux ratios have a logarithmic 
standard deviation corresponding to a flux ratio of 2.5, and 
a flux ratio of more than 10 corresponds to a 2.5er level of 
variability. We label this class “continuum variability” and note 
that the variability arises from a wide range of causes, including 
most AGN variability, measurement errors, and differences in 
the spectral response of the two instruments. 

A different, more dramatic, type of variability was also 
present in the survey. At high flux values (toward the right 
in Figure 2), there is a clear separation between the sources 
which were observed in the RASS, and those which were 
not observed in the RASS. The ratios between the XMMSL1 
flux and the RASS upper limit for sources not detected in 
RASS are plotted as red triangles in the figure. This second 
population (state-change objects) represents sources which were 
in some distinctly different state at the time of the RASS 
and XMMSL1 observations. Known examples in this class 
include X-ray binaries in different states, flaring stars, tidal 
disruption events (TDEs), and classical novae. This class could 
be defined as sources with a flux and flux ratio which places 
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Table 1 

The List of Transient Objects with Galaxy Associations Found in our Search 


Name 

Flux 

Ratio 


Luminosity 

Notes 

Date 

XMMSLl... 

(cgs) 


(Mpc) 

(ergs -1 ) 



J202320.7— 670021 

2.3 x 10 -11 

46 

67 

1.2 x 10 43 

Not AGN C 

2009 Oct 05 

J084837. 9+193527 

6.4 x 10~ 12 

23 

283 

6.1 x 10 43 

Emission line galaxy d 

2010 Oct 17 

J 182609.9+545005 

4.5 x 10~ 12 

22 

650 

2.2 x 10 44 

Not AGN 

2005 Nov 01 

J152408. 6+705533 

3.3 x 10~ 12 

15 

256 

2.6 x 10 43 

Spiral galaxy 

2006 Jan 24, 2007 Nov 01 

J202554.8-5 11629 

9.1 x 10- 12 

12 

? 

? 


2010 Apr 16 

J 1 3 1 95 1 .9+225957 

3.7 x 10~ 12 

10 

99 

4.3 x 10 42 

Not AGN d 

2005 Jul 15 

J1 11527.3+180638 

7.1 x 10- 12 

30 

12 

1.2 x 10 41 

Esquej. 

2003 Nov 22 

J155631. 5+632540 

3.8 x 10~ 12 

21 

? 

? 

Bad match 

2010 Jun 23 

J020303. 1—074154 

3.1 x 10~ 12 

12 

272 

2.7 x 10 43 

Esquej. 

2004 Jan 14 

J 1 70543.0+850523 

4.5 x 10~ 12 

12 

? 

? 

Cluster 

2008 Aug 31 

J013727.9- 195605 

1.2 x 10 -11 

11 

1320 

2.1 x 10 45 

Cluster 

2007 Dec 27 

J104745. 6— 375932 

3.6 x 10~ 12 

10 

335 

4.8 x 10 43 

AGN C 

2003 Dec 14 


Notes. The listed fluxes are those reported in XMMSL1 catalog (Saxton et al. 2008). The column labeled “Ratio” shows the ratio between the observed XMMSL1 
flux and a 2a upper limit based on the corresponding RASS data. Each distance (Di) and inferred luminosity is based on the overlap of the 2a XMMSL1 position 
with an optically identified galaxy. The note “Not AGN” signals that an available spectra shows no evidence for AGN emission. The objects with the note Esquej were 
reported as tidal disruption event candidates by Esquej et al. (2007). Bad match denotes that the matched galaxy was separated from the XMMSL1 source by more 
than the 2a position uncertainty. The six objects above the horizontal line are further discussed in Section 3. 
a The luminosity distance of the associated galaxy. 
b Based on inspection of 6dF spectrum. 
c Based on inspection of SDSS spectrum. 
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Figure 3. Histogram of the logarithm of the flux ratio for objects detected in 
both XMMSL1 and RASS. The standard deviation of the distribution is 0.40, 
corresponding to a flux ratio of 2.51. By this measure, then, a threshold on 
the flux ratio of 10 requires an object’s variability to be beyond 2.5a in the 
distribution. 

(A color version of this figure is available in the online journal.) 

them outside the distribution of continuum variability shown in 
Figure 3. Many of the sources with RASS non-detections (red 
triangles) seem to belong to this class. Though the separation 
between the RASS detections and non-detections in Figure 2 
disappears at low flux values, the vertical placement of the non- 
detections is determined primarily by the limiting flux of the 
RASS observation. Therefore, any of the non-detections would 
potentially sit higher in this plane with data from a deeper 
observation. Given that some of the bright XMMSL1 sources 
were seen to vary by more than a factor of 100, it would be 
surprising if this was not also the case for some of the dimmer 
sources as well. 

Some of the state-change objects observed in XMMSL1 were 
also studied by Starling et al. (201 1). In their study, the authors 


selected 97 sources from v 1 .4 of the XMM-Newton Slew Survey 
with no identified optical or RASS counterpart. They then 
collected pointed observations for 94 of these with Swift /X RT 
in an attempt to classify them. In their paper, Starling et al. 
(201 1) report that 71% of the targets were not seen in the Swift 
observations, implying they had faded by at least a factor of 
10. The isotropic distribution and lack of an optical counterpart 
of these mysterious sources led the authors to conclude that 
they represented a primarily extragalactic population. None 
of the sources in the Starling study appeared in our final 
candidate list, as the former campaign required sources to have 
no optical counterpart, where we required a match to a known 
optical galaxy. However, if the bulk of the transients studied 
by Starling et al. were extragalactic, as they concluded, then 
the objects listed in Table 1 may represent low redshift objects 
from the same population. In particular, the study found 47 
objects detected in the XMM-Newton soft band, but not detected 
with Swift. 

We wished to use the variability in this data set to characterize 
the chance of a spurious detection in coincidence with a trigger 
from the future advanced LIGO/ Virgo network. We applied 
the selection criteria described in Section 2, and the statistical 
results of our search are presented in Figure 4. Each curve in 
the figure shows the density of sources in the set with fluxes 
above the value on the X-axis. The top curve (black, dotted) 
is the log /V-log S curve of the XMMSL1 catalog, assuming 
32,800 deg 2 of coverage. The bottom curve (blue, solid) is the 
result of our search criteria, which resulted in 1 2 transient objects 
spatially coincident with a known galaxy (see Table 1, yellow 
circles in Figure 2). In requiring the galaxy association, we 
included associations with both galaxies and galaxy clusters, but 
rejected associations labeled in the XMMSL1 as AGN (Seyfert, 
BL Lac, etc.). 

Most of the matched host galaxies with available redshifts 
were located at luminosity distances of less than 350 Mpc. 
This is most likely because available catalogs of galaxies were 
dominated by relatively shallow observations from large area 
surveys. Kasliwal (2011) found that current galaxy databases 
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Figure 4. The statistical result of our search for low-redshift transients. The top 
curve (black) shows the log/V-logS" plot for the XMMSL1 catalog, containing 
1411 objects above our flux threshold. The second curve (red) shows the 
distribution of the 97 transients, defined as at least 10 times brighter in XMMSL1 
than in RASS. The bottom curve includes only transient objects that are spatially 
coincident with a known galaxy, after rejecting previously identified AGNs. 

(A color version of this figure is available in the online journal.) 

(NED, 10 HyperLeda, 11 etc.) were 50% complete to 200 Mpc, 
with a steep downward trend in completeness as a function 
of distance. This distance scale was well aligned with the 
reach of Advanced LIGO and Virgo, which will detect NS-NS 
mergers primarily between 100 and 400 Mpc (Nissanke et al. 
2013). From an operational perspective, then, future searches for 
X-ray counterparts to LIGO/ Virgo transients may not need to 
obtain a redshift estimate for each possible host galaxy in the 
field; any galaxy bright enough to be cataloged by a large 
area survey without an AGN signature is likely to be within 
the Advanced LIGO horizon. For our interpretation, we have 
taken this approach. This would tend to skew our estimates 
of unrelated counterparts to be artificially high, because using 
the estimated distance to the gravitational wave source can be a 
powerful tool for rejecting potential host galaxies at inconsistent 
redshifts (Nissanke et al. 2013). 

To interpret Figure 4, we note that typical exposures in 
the XMM-Newton Slew Survey were 5-20 s, so we assumed 
that the transients that passed our selection criteria were of 
longer duration than the XMM-Newton exposures. Though the 
Slew Survey was not uniform in any sense, the coverage area 
was largely random, and so we assumed that the density of 
sources which passed our cuts was not strongly biased by 
target selection. The survey did include some repeat visits, and 
only covered 20,900 unique square degrees. However, the time 
between repeated visits was typically > 1 yr, which was more 
than the expected timescale for fading of X-ray counterparts to 
neutron star mergers. For this reason, we interpreted our results 
based on 32,800 deg 2 of coverage, noting that this choice leads 
to at most a 60% systematic in our results. So, for example, 
the 12 sources above our flux threshold may be interpreted as 
4 x 10” 4 transients per square degree on timescales shorter than 
the ~15 yr between RASS and XMMSL1. 

The positional uncertainties associated with a trigger from 
the LIGO/ Virgo network are expected to vary a great deal 


10 http://ned.ipac.caltech.edu/ 

11 http://leda.univ-lyonl.fr/ 


depending on signal-to-noise ratio and other factors; however, 
studies typically quote numbers between 20 and 200 deg 2 for 
the position uncertainty of a low signal-to-noise ratio trigger 
with three sensitive gravitational wave detectors operational 
(Abadie et al. 2012; Fairhurst 2009; Klimenko et al. 2011). 
Searching such a large area for an X-ray counterpart will require 
a relatively high flux limit to the search, since the instrument 
will have to either be very wide field or will have to use short 
exposures for many tiles. For example, in principle Swift /X RT 
could utilize ~100 s exposures to tile as much as 35 deg 2 in 
a day, though only to a depth of 6 x 10” 12 erg cirT 2 s” 1 
(Kanner et al. 2012). While practical issues concerned with 
many repointings of the instrument may make this difficult in 
practice, Figure 4 shows that there would be less than a 1% 
chance of finding a transient in this search area by chance. A 
more natural scenario would be the application of a very wide- 
held, focusing instrument such as the proposed ISS-Lobster. 
The proposed instrument would image a 400 deg 2 held of 
view to a depth of around 10” 11 erg s” 1 cm” 2 in a 20 minute 
exposure. Applying these numbers, this observation would have 
a 3% chance of imaging an unrelated transient coincident with a 
known host galaxy, or less than a 1% chance if we imagine only 
using the fraction of the held that overlaps the LIGO/ Virgo 
errorbox. Similar considerations would apply to other wide 
held-of-view instruments, including MAXI (Matsuoka et al. 
2009) or the proposed A-STAR (Osborne et al. 2013). 

Of course, the details of how a future search is carried out 
would have a strong inhuence on these numbers. Choosing 
to dehne a transient as an object that brightens by at least a 
factor of 10 is somewhat arbitrary, though for our data set, 
the naturally distinct classes of continuum variability and state- 
change variability seen in Figure 2 appears to somewhat justify 
this choice. Another important factor is the completeness of 
the galaxy catalog, which is likely to evolve rapidly due to 
efforts by several large area surveys (Metzger et al. 2013). On 
the other hand, at the order-of-magnitude level, these results 
seem to be robust. The associations of these sources with the 
host galaxies do not appear to be spurious (See Section 4.5), 
so adjusting the match radius used to associate host galaxies 
and sources will make little impact. Similarly, increasing the 
searched area around the host galaxy due to “kicks” in the 
binary (Fryer & Kalogera 1997) will not change these statistics 
substantially. A variety of models, with some validation from 
studies of short GRB host galaxies, have concluded that the 
majority of NS-NS mergers occur within 10 or a few tens of 
kiloparsecs from the centers of their host galaxies (Berger 2010; 
Brandt & Podsiadlowski 1995; Bloom et al. 1999; Fryer et al. 
1999; Belczynski et al. 2006). At 200 Mpc, the 30" search radius 
used for galaxy association corresponds to an offset from the 
host galaxy of 30 kpc. Berger (2010) showed that for a range of 
models, this search radius would include 70%-90% of NS-NS 
mergers. On the other hand, models do exist with more extreme 
kick velocities, leading to mergers that occur up to a megaparsec 
away from the host galaxy (Kelley et al. 2010). To accommodate 
these models would mean using a somewhat larger search radius, 
and so the requirement of a galaxy association may become less 
useful. Additional surveys could increase the number of known 
galaxies, and so increase false associations by a factor of ~2 
or more. However, eliminating possible hosts with redshifts 
inconsistent with the GW data could limit this effect (Nissanke 
et al. 2013). 

Finally, we note that the separation of soft X-ray sources 
into state-change transients and continuum variability appears to 


4 


The Astrophysical Journal, 774:63 (9pp), 2013 September 1 


Kanner et al. 



Figure 5. Optical spectrum of the host galaxy matched to XMMSL1 
J104745. 6— 375932 obtained by the 6dF survey, plotted in the restframe us- 
ing a redshift of 0.075. The strong, broad hydrogen lines and the blue spectrum 
are characteristic of quasars. 

(A color version of this figure is available in the online journal.) 



Figure 6. Optical spectrum of the host galaxy matched to XMMSL1 
J202320.7— 670021, obtained by the 6dF survey, and plotted in the restframe 
using a redshift of z = 0.016. The galaxy is dominated by stellar absorption 
features. 

(A color version of this figure is available in the online journal.) 


occur naturally. For these reasons, we expect that the finding that 
around 10% of bright, soft X-ray sources demonstrate a state- 
change brightening over long timescales, and that around 10% of 
these can be associated with galaxies within 200 Mpc, will prove 
true for future searches with perhaps a factor of a few uncertainty 
in both cases. For searches for X-ray counterparts to advanced 
LIGO/ Virgo triggers, this criteria represents a two order of 
magnitude reduction in the background rate, as compared with 
the density of X-ray sources on the sky which has been used 
to estimate the density of spurious associations in past work 
(Evans et al. 2012). 

3.2. Further Selection 

After identifying our list of 12 candidates, we sought to char- 
acterize them as far as possible. We inspected optical images 
of the host galaxies and searched for additional data using the 
FIEASARC database. We also inspected publicly available and 
newly obtained optical spectra to check for AGN signatures. 
Two of the objects were rediscoveries of previously published 
candidate TDEs (XMMSL 1 J020303 .1-074154 and XMMSL 1 
J1 11527.3+180638; Esquej et al. 2007). One object showed 
broad AGN emission lines in a 6dF 12 archived spectrum, and 
was seen in XRT data to have an X-ray spectrum consistent with 
an AGN (XMMSL1 J104745.6-375932; see Figure 5). One ob- 
ject’s (XMMSL 1 J 1 5563 1 .5+632540) 2 o position circle did not 
include the matched host galaxy. Two of the sources (XMMSL 1 
J013727.9- 195605 and XMMSL1 J170543.0+850523) were 
matched to galaxy clusters with ambiguous galaxy associations. 

This left six objects which exhibited some interesting prop- 
erties. Of the five where we obtained a redshift (either through 
observations or the literature), most corresponded to luminosity 
distances less than 300 Mpc in a standard cosmology; the fur- 
thest was placed around 650 Mpc. We examined optical spectra 
for five of the host galaxies, and none of these showed broad 
emission lines, though two (XMMSL1 J084837. 9+193527 
and J152408. 6+705533) showed narrow emission lines (see 
Figures 6-10). They were all soft sources, and only one has 
a measured hardness ratio presented in XMMSL 1 (XMMSL 1 


12 http://www.aao.gov.au/6dFGS 



Figure 7. Optical spectrum of the host galaxy matched to XMMSL 1 
J084837.9+ 193527, obtained by the SDSS, plotted in the rest frame using a 
redshift of 0.064. A test using diagnostic narrow line ratios shows this to be a 
star forming galaxy, as described in the text. 

(A color version of this figure is available in the online journal.) 

J 182609.9+545005: -0.46). We searched HEASARC for ob- 
servations with Chandra, Swift, and XMM-Newton, but found 
no observations containing any of the six remaining sources. 
These six objects were generally extragalactic, more powerful 
than4x 10 42 erg s -1 , close (Di < 300 Mpc), variable by at least 
lOx, and lacking evidence for AGN activity. This makes them 
difficult to characterize, as well as energetic and variable. In 
the next section, we discuss possible characterizations for these 
sources, and attempt to address the likelihood of each. 

4. POSSIBLE CHARACTERIZATIONS 
4.1. Tidal Disruption Events 

TDEs are thought to occur at a rate of ~ 1 0 5 yr -1 Mpc -3 
and display initial luminosities up to ~10 45 erg s during the 
first few days or weeks followed by a characteristic dimming 
oc f _5 /3 (Bade et al. 1996; Wang & Merritt 2004). Under these 
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Figure 8. Optical spectrum of the host galaxy matched to XMMSL1 
J 1 82609.9+545005, obtained with the Keck/DEIMOS on 2013 May 3, plot- 
ted in the rest frame with a redshift of 0.14. The absence of emission lines and 
red color suggest this is a red galaxy without an active central region. 

(A color version of this figure is available in the online journal.) 



Figure 10. Optical spectrum of the host galaxy matched to XMMSL1 
J1 3 195 1.9+225957, obtained by the SDSS, plotted in the rest frame with a 
redshift of 0.023. The spectrum appears to be dominated by stellar absorption 
features. 

(A color version of this figure is available in the online journal.) 



Figure 9. Optical spectrum of the host galaxy matched to XMMSL1 
J1 52408.6+705533, obtained on 2013 April 15 with the 200 inch Hale Tele- 
scope at Palomar Observatory, plotted in the rest frame with a redshift of 0.059. 
The narrow line features and WISE colors suggest this is a spiral galaxy, and so 
is unlikely to host an active nucleus. 

(A color version of this figure is available in the online journal.) 

assumptions, at 200 Mpc, such an event would be visible for a 
few months to years above a flux of 3 x 10~ 12 erg s 1 cuT 2 . A 
number of TDE candidates have been discovered through their 
X-ray emission (Komossa & Bade 1999; Grupe et al. 1999; 
Komossa & Greiner 1999; Greiner et al. 2000b; Maksym et al. 
2010; Lin et al. 201 1). In addition, there have been discoveries 
made using ultraviolet surveys, such as those discussed in Gezari 
et al. (2006, 2012), and with optical surveys (van Velzen et al. 
2011; Drake et al. 201 1). There have also been two well-studied 
TDEs discovered with the Swift satellite (e.g., Bloom et al. 2011; 
Levan et al. 2011; Burrows et al. 2011; Cenko et al. 2012). 

In fact, some TDEs have already been discovered using data 
from the XMM-Newton Slew Survey. Last year, Saxton et al. 
(2012) reported on the discovery of an X-ray flare in SDSS 
J 120136. 02+300305. 5 found in data from the XMM-Newton 
Slew Survey, and followed up with pointed observations of 


Table 2 

Galaxy Matches from the XMMSL1 Catalog for Six Transients 
Identified by our Search 


Name 

XMMSL1. . . 

Galaxy Match 

Offset 

(arcsec) 

J202320.7— 670021 

6dFGS gJ202322.7— 670046 

27 ± 13 

J084837.9+193527 

SDSS J084838. 57+193528.9 

10 + 22 

J 182609.9+545005 

2MASX J1 8261094+5450052 

7 + 5 

J152408.6+705533 

MRK 1097 

1 ±3 

J202554.8— 511629 

2MASX J20255579-51 16276 

9 + 19 

J13 195 1 .9+225957 

NGC 5092 

4 + 2 


Notes. The offset shows the angular distance between the galaxy and the source 
position, along with the uncertainty quoted in the catalog. 


Swift/XRT. Their program is designed to find flaring X-ray 
events soon after they begin to enable prompt follow-up. The 
flux of the corresponding observation listed in the XMMSL1 
was below our flux threshold, so this event was not rediscovered 
by our search. Esquej et al. (2007, 2008) found two TDEs in the 
XMM-Newton Slew Survey, including one TDE above our flux 
threshold (see Table 1) , via a search with similar selection 
criteria, but using the first release (vl.l) of the XMMSL1 
catalog, covering 6300 deg 2 . Our search used the latest release 
(vl.5) of XMMSL1, covering 32,800 deg 2 (20,900 deg 2 of 
unique area). This suggests that we should expect around four 
new TDEs in our data set. Given the lack of nuclear activity in 
some of the host galaxies, the high luminosity of the sources, 
and the positions which are consistent with the center of the host 
galaxies, we expect this model will account for at least some 
of our candidates. For example, XMMSL1 J13 195 1 .9+225957 
was matched to NGC 5092 with an offset from the center of 
the galaxy of only 4", compared with the match radius of 30", 
or the average XMMSL1 uncertainty of 8" (see Table 2). An 
SDSS spectrum of this galaxy, shown in Figure 10, reveals only 
absorption features, making an AGN association unlikely. These 
features were consistent with a TDE description. The galaxy 
SDSS J084838. 57+193528. 9 showed some narrow emission 
line features. To characterize it, we used the best fit line profiles 
provided by the SDSS Science Archive Server. We found the 
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[N ii] A6583/Ha ratio to be 0.32, and the [O m] A5007/H/3 ratio 
to be 0.67, so that this galaxy appeared to have a star-forming 
region, but no nuclear activity (Veilleux & Osterbrock 1 987). For 
this reason, the observed variability in this galaxy also appeared 
consistent with a tidal disruption model. We obtained a spectrum 
of the host galaxy matched to XMMSL1 J 182609.9+545005 on 
2013 May 3, using the Keck/DEIMOS (Figure 8). The spectrum 
was dominated by absorption features, and showed no evidence 
for an active central region. The redshift we obtained (z = 0.14) 
places this galaxy outside the planned Advanced LIGO horizon. 

Finally, one object that passed our cuts, XMMSL1 
J152408. 6+705533, was observed twice in the XMMSL1, in 
2006 January and 2007 November. In the 23 months between the 
observations, the source faded by a factor of three. Fitting a TDE 
light curve to these two points resulted in an event with starting 
time in 2004 March, and a luminosity of 5.7x 10 43 erg s” 1 1 yr 
after the start time. The implied energies were roughly consis- 
tent with previously observed TDEs (Esquej et al. 2008). We 
obtained a spectrum of the host galaxy, MRK 1097, with the 
200 inch telescope at Palomar observatory, and found primarily 
narrow line emission features (see Figure 9). The WISE col- 
ors for this galaxy ([W l — W2] = 0.3 mag and [W 2 — W 3] = 
1.8 mag) (Wright et al. 2010) and narrow line features suggested 
this is a spiral galaxy, and so is unlikely to host an AGN. 

4.2. GRB Afterglows 

Short GRBs have been observed with a rate density of 
5-13 Gpc” 3 yr” 1 , implying a rate of around one per year within 
300 Mpc (Nakar 2007; Coward et al. 2012). Long GRBs are 
observed with a rate density of 0.5 Gpc -3 yr” 1 (Nakar 2007), 
or one event every 20 yr within 300 Mpc. GRBs are known to 
display bright afterglows, typically observable in soft X-rays 
for a few hours up to a few days after the burst. However, most 
GRBs are observed much further away than 300 Mpc, so a 
burst at the distance of our objects would have an afterglow 
that would appear brighter, and would potentially be observable 
longer. Taking an optimistic but plausible scenario, a short 
GRB afterglow at the distances our objects, showing power-law 
dimming with a temporal index of 1.2, would display the flux 
levels observed with XMM-Newton 10-100 days after the burst. 
A model that includes a jet break would show a faster fading, 
with a temporal index closer to 2.0 (Racusin et al. 2009). This 
means that, even in the optimistic case, we might only expect 
one GRB per year within our search volume, and that a GRB’s 
afterglow would only be observable for one to three months. 
Based on these numbers, we believe there is less than a ~ 1 0% 
chance that our candidate list, after the requirement of a galaxy 
association, includes one or more GRB afterglows. 

On the other hand, there are related classes of transients, 
both observed and theoretical, that are thought to be more 
common in the local universe. Low-luminosity GRBs have 
been observed with a local rate density much higher than 
the rate of cosmologically observed bright GRBs (Soderberg 
et al. 2006; Cobb et al. 2006; Chapman et al. 2007). These 
events can have X-ray band afterglows that are less luminous 
than cosmological GRB afterglows, but otherwise with similar 
properties (Soderberg et al. 2006). So-called orphan afterglows 
and failed GRBs may also be more common in the local universe 
than the more commonly observed cosmological GRBs (Rhoads 
2003; Huang et al. 2002), an idea that may be supported by a 
recent observation (Cenko et al. 2013). It is difficult to rule out 
the possibility that such an event is in our sample. 


4.3. Ultra-luminous X-ray Source 

Ultra-luminous X-ray sources (ULXs) have been ob- 
served in several nearby galaxies, with X-ray luminosities of 
10 38 — 10 41 erg s” 1 . These objects are known to exhibit short 
time variability, and have been observed in both high and low 
energy states (Winter et al. 2006). However, the luminosities 
of our sources exceed the range of known ULXs, so, these ob- 
jects cannot be naturally characterized as ULXs, or at best, they 
would represent extreme examples of the class. 

4.4. AGN 

We inspected optical spectra for five of our six objects, none 
of which showed evidence for an active central region. Two 
of the spectra showed no strong emission lines, while two 
host galaxies (SDSS J084838.57+193528.9 and MRK 1097) 
showed star formation lines. It is possible that one or both of the 
galaxies for which we do not have optical spectra will turn out 
to be an AGN. Moreover, there are known cases of galaxies that 
seem consistent with an AGN model when observed in X-rays, 
but do not show evidence of nuclear activity in their optical 
spectra (Jackson et al. 2012). The details of the mechanism that 
hides the active region is still being disputed. While AGNs are 
known to exhibit variability, both on long and short timescales, 
the majority of this variability is low amplitude. For example, 
Saxton et al. (2011) found that, in a sample of over 1000 AGNs 
observed with both the XMM-Newton Slew Survey and RASS, 
only 5% varied by more than a factor of 10. Given the relative 
rarity of large amplitude variability in AGNs, and the current 
difficulty in describing optically quiescent but X-ray luminous 
galaxies, a low redshift, variable, “hidden” AGN might be an 
interesting source for future study. 

4.5. Spurious Detections 

Given that we have selected objects that appeared in 
XMMSL1, but not in RASS, one has to consider the possibility 
that these are spurious detections, or perhaps that they are real 
sources, but the galactic associations are incorrect. We note that 
XMMSL1 is estimated to contain less than 1% false sources in 
the soft band (Saxton et al. 2008), though the stronger argument 
is that the associations with host galaxies are unlikely to all be 
false. To make a conservative estimate of the chance of false 
coincidences, we assumed a 30" match radius, and imagined 
spreading the 94 objects in our sample with flux ratio greater 
than 10 randomly across the sky, so that there were 1.7 x 10” 3 
transient objects deg” 2 . We found that the odds of finding the 
observed number of coincidences by chance were very small. 
For example, within 200 Mpc there are ~ 150,000 known galax- 
ies (3.6 galaxies deg” 2 ) (Kasliwal 2011). So, within 200 Mpc, 
our search had a 7% chance of finding even a single coincidence 
by chance, where we found four. Moreover, all but one of the 
sources matched their host galaxies within 15" despite a 30" 
search radius, and none of the transients were spatially incon- 
sistent with having originated from within their hosts’ angular 
extent (see Table 2). For these reasons, most or all of the claimed 
galaxy associations are very likely real. It is worth noting, how- 
ever, that one object has a slightly larger offset from the host 
galaxy, namely J202320. 7— 670021. The matched host galaxy 
is at a distance of only 67 Mpc, and the angular extent of the 
galaxy overlaps the ler error circle from XMMSL1. So, while 
it seems likely that this association is real, in this case we are 
unable to rule out a spurious association. 
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5. CONCLUSIONS 

In this work, we took a census of X-ray transient objects 
in the low redshift universe, motivated by future observa- 
tions of compact object mergers with Advanced LIGO and 
Advanced Virgo. A wide-field, soft X-ray monitor such as 
the proposed ISS-Lobster will be able to seek counterparts 
to LIGO/ Virgo events, but high-confidence identifications will 
demand high quality variability studies. We performed a sys- 
tematic search for low-redshift, extragalactic transients in the 
XMM-Newton slew survey, covering 32,800 deg 2 , above fluxes 
of 3 x 10 -12 erg cm -2 s -1 in the 0.2-2 keV band. We com- 
pared observations taken many years apart in the RASS and 
the XMMSL1, and found that over these long timescales, vari- 
ations in flux divided sources into two categories: continuum 
variability, which can be described as a log-normal distribution 
in flux variability with a width of around 3, and state-change 
variability, corresponding to objects which show more dramatic 
changes in flux. State-change sources at low redshift may be 
confusion sources for future searches for counterparts to events 
measured with Advanced LIGO and Advanced Virgo, so we 
sought to characterize their density on the sky. We found that 
transient sources represented around 10% of all objects in a flux 
limited survey, and that of these, around 10% could be associ- 
ated with known optical galaxies. For searches for LIGO /Virgo 
counterparts using wide-field X-ray imagers capable of observ- 
ing hundreds of square degrees, we should expect around one 
false coincidence for every 10,000 deg 2 searched to a flux limit 
of 10 _n erg s -1 cm -2 . 

We found 12 objects meeting our search criteria, most of 
which were located within 350 Mpc. Of these, we identified six 
with clear galaxy identifications that were difficult to character- 
ize and that have not been previously studied. They were highly 
luminous, highly variable, and lacked classical AGN optical 
signatures in the five cases with available spectra. Four of the 
sources (XMMSL1 J13 195 1 .9+225957, J084837.9+193537, 
J 182609.9+545005, and J 152408.6+705533) met all of the cri- 
teria of Esquej et al. (2007) for identifying candidate TDEs. 
One other seemed consistent with a TDE description, but may 
prove to be a variable AGN, and the sixth source had a measured 
position only marginally consistent with the center of the appar- 
ent host galaxy. It is possible that, whatever the nature of these 
sources, they represent some of the closest members of the class 
of unexplained, extragalactic transients identified by Starling 
et al. (2011), who focused on transient objects not associated 
with known galaxies. This study represents the first attempt to 
characterize soft X-ray confusion sources for counterparts to 
Advanced LIGO/ Virgo merger events, and in the process, has 
revealed a handful of unusual objects that are both powerful and 
dynamic. 
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